Non-Abelian topological phases in an extended Hubbard model 
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We describe four closely related Hubbard-like models (models A,B,C and D) of particles that can hop on a 
2D Kagome lattice interacting via Coulomb repulsion. The particles can be either bosons (models A and B) or 
(spinless) fermions (models C and D). Models A and C also include a ring exchange term. In all four cases 
we solve equations in the model parameters to arrive at an exactly soluble point whose ground state manifold is 
the extensively degenerate "rf-isotopy space" V < d < 2. Near the "special" values, d = 2cos7t/fc + 2, 
should collapse to a stable topological phase with anyonic excitations closely related to SU(2) Chern-Simons 
theory at level k. We mention simplified models A~ and C which may also lead to these topological phases. 



I. INTRODUCTION 

Since the discovery of the fractional quantum Hall ef- 
fect in 1982i, topological phases of electrons have been 
a subject of great interest. Many abelian topologi- 
cal phases have been discovered in the context of the 
quantum Hall regime 2 . More recently, high-temperature 
superconductivity 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 and other com- 
plex materials have provided the impetus for further theoreti- 
cal studies of and experimental searches for abelian topolog- 
ical phases. The types of microscopic models admitting such 
phases is now better understood 17 - 18 - 19 - 20 . 

Much less is known about non-abelian topological phases. 
They are reputed to be obscure and complicated, and there 
has been little experimental motivation to consider non- 
abelian topological phases, apart from some tantalizing hints 
that the quantum Hall plateau observed at v = 5/2 might 
be non-abeliar. 21 * 2 ^ 2 ^ 4 * 2 ^ 2 ^ 2 ^. However, non-abelian topo- 
logical states would be an attractive milieu for quantum 
computation 28,29 . Furthermore, a better understanding of non- 
abelian topological phases would help determine if they have 
already been unwittingly observed in nature. It is, therefore, 
encouraging that recent progress in understanding the Hamil- 
tonian formulation of a class of topological field theories 
has reduced computations in these theories to combinatorial- 
pictorial manipulations^^!. These theories do not break par- 
ity or time-reversal symmetry, which may have practical ad- 
vantages. A further virtue of this formulation of these theo- 
ries is that it exposes a strategy for constructing microscopic 
physical models which admit the corresponding phases; since 
Hilbert space is reduced to a set of pictorial rules, the mod- 
els should impose these rules as energetically-favorable con- 
ditions satisfied by the ground state. In this paper, we show 
how this approach can be implemented. 

The microscopic models which we construct are 'quasi- 
realistic'. Their precise form is not likely to be realized in 
any material. However, they are soluble, so we know that 
they do, indeed, support the non-abelian topological states of 
matter which interest us, and one can imagine a real material 
whose Hamiltonian can be viewed as a small perturbation of 
one of the Hamiltonians of this paper. It may also be possible 
to design a quantum dot or Josephson junction array with such 
a Hamiltonian. 



The non-abelian topological phases which arise in this pa- 
per are related to the doubled SU(2)^ Chern-Simons theories 
described in Refs. l3ll32l In appendixlAlsome further informa- 
tion about these theories, called (DK^, is giver. 48 . For the mo- 
ment, we note that (DK^ has ground state degeneracy (k+ l) 2 
on the torus T 2 and should be viewed as a natural family con- 
taining the topological (deconfined) phase of Z2 gauge theory 
as its initial element, k = 1. For k > 2 the excitations are 
non-abelian. For k = 3 and k > 5 the excitations are computa- 
tionally universal 3 ^ 4 .. In Appendix [B] some general purpose 
formulae for the perturbation theory of Hubbard models are 
derived. 

The basic strategy of the present paper is to construct a 
Hamiltonian which enforces a subset - c/-isotopy 31 - 32 - of the 
defining conditions of the Hilbert space of these theories. The 
resulting low-energy subspace will be called Vd- It has been 
argued 31 - 32 that for special values of the parameter d, namely 
d = 2cos(n/k + 2), Vd will collapse to the stable topological 
phase DK^. 



II. THE BASIC MODEL 

The models described in this paper, both bosonic and 
fermionic, are defined on the Kagome lattice. The basic lattice 
is shown in Fig. [2 The sites of the lattice are not completely 
equivalent, in particular we choose two special sublattices - 
^ (red) and Q (green) whose significance shall became clear 
later on. In some cases we will introduce "defects" in this 
sublattice arrangement. 

Our generic Hamiltonian is given by: 



(ij)eO {i,j)ex40 

Y J t i j(c]cj+c]c,)+ Ring. (1) 

(U) 



Here = cjc; is the occupation number on site i, ^ is the cor- 
responding chemical potential. Uo is the usual onsite Hubbard 
energy Uo (clearly superfluous for spinless fermions). U is 
a (positive) Coulomb penalty for having two particles on the 
same hexagon while Vy represent a penalty for two particles 
occupying the opposite corners of "bow ties" (in other words, 
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FIG. 1: Kagome lattice with the special sublattices (depicted 
red) and g (green). 



being next-nearest neighbours on one of the straight lines). 
We allow for the possibility of inhomogeneity so not all Vy 
are assumed equal. Specifically, define v c uh = Vy where a is 
the color of site (z), b is the colour of (/'), and c is the colour 
of the site between them. In the lattice in Fig.^ we have, pos- 
sibly distinct, v^ fe , v* fo , vf , v h rb , and v b rg , where r e g e Q, 
and b € S = i*C\(X U f, ; - is the usual nearest-neighbour 
tunnelling amplitude which is also assumed to depend only on 

; f , where c now refers to 



sublattice themselves form one such dimer covering, a so- 
called "staggered configuration". This particular covering is 
special: it contains no "flippable plaquettes", or rhombi with 
two opposing sides occupied by dimers (see Fig.|2}. 



the colour of the environment: t, 



ij 

the colour of the third site in a triangle. Finally, "Ring" is a 
ring exchange term - an additional kinetic energy term which 
we add to the Hamiltonian on an ad hoc basis to allow cor- 
related multi-particle hops which "shift" particles along some 
closed paths (see more on this term below). Ring exchange 
terms can be justified semiclassically, but are not generally in- 
cluded as bare terms in the Hubbard model 49 . In models B and 
D we remove the Ring term at the expense of complicating 
the sublattice (and thus the interactions) with the addition 
of several new colors. 

The onsite Hubbard energy Uq is considered to be the 
biggest energy in the problem, and we shall set it to infin- 
ity, thereby restricting our attention to the low-energy mani- 
fold with sites either unoccupied or singly-occupied. The rest 
of the energies satisfy the following relations: U ^> fy, Vy,;U;; 
we shall be more specific about relations between various fy's, 
Vi/s and /j,'s later. 

We are particularly interested in the 1/6-filled case (i.e. 
N p = = N/6, where N is the number of sites in the lat- 
tice). The lowest-energy band then consists of configurations 
in which there is exactly one particle per hexagon, hence all 
U -terms are set to zero. (One example of such a configura- 
tion is given by particles occupying all sites of sublattice ^..) 
These states are easier to visualise if we consider a triangular 
lattice 1 whose sites coincide with the centers of hexagons of 
3(\ (ac is a surrounding lattice for 1 .) Then a particle on 
is represented by a dimer on T connecting the centers of two 
adjacent hexagons of The condition of one particle per 
hexagon translates into the requirement that no dimers share a 
site. In the 1/6-filled case this low-energy manifold coincides 
with the set of all dimer coverings (perfect matchings) of 1 . 
The "red" bonds of T (the ones corresponding to the sites of 




FIG. 2: Triangular lattice T obtained from ^ by connecting the 
centers of adjacent hexagons. The bonds corresponding to the special 
sublattices %, and Q are shown in red and green correspondingly. See 
text for more on bond color-coding. Triangles with one red side are 
shaded as guide to the eyes: these triangles will be essential for the 
dynamics of our models. 

So henceforth particles live on bonds of the triangular lat- 
tice (Fig. |3 and are represented as dimers^. In particular, a 
particle hop corresponds to a dimer "pivoting" by 60° around 
one of its endpoints, Vu — v c ab is now a potential energy of 
two parallel dimers on two opposite sides of a rhombus (with 
c being the color of its short diagonal). 

From the relation to dimer coverings exhibited above, it is 
clear that our model is in the same family as the quantum 
dimer model 35 , which has recently been shown to have an 
abelian topological phase on the triangular lattice 17 which, in 
the notation of this paper, is <DK\. Here, we show how other 
values of k can be obtained. 

The goal now is to derive the effective Hamiltonian acting 
on this low-energy manifold represented by all possible dimer 
coverings of T . Our analysis is perturbative in t/U =: e. 

The initial, unperturbed ground state manifold for Uq = °°, 
U large and positive, all f i; - , V !; - = and all /j, equal is spanned 
by the dimerizations © of the triangular lattice 1 . As we grad- 
ually turn on the f's, v's, and Ring terms, we shall see what 
equations they should satisfy so that the effective Hamiltonian 
on © has the desired space Vd as its ground state manifold 
(GSM). 

Since a single tunnelling event in © always leads to dimer 
"collisions" (two dimers sharing an endpoint) with energy 
penalty U, the lowest order at which the tunnelling processes 
contribute to the effective low-energy Hamiltonian is 2. At 
this order, the tunnelling term leads to two-dimer "plaquette 
flips" as well as renormalisation of bare onsite potentials fjfs 
due to dimers pivoting out of their positions and back. We 
always recompute bare potentials /j,'s to maintain equality up 
to errors (e 3 ) among the renormalized /Vs. This freedom to 
engineer the chemical potential landscape to balance kinetic 
energy is essential to finding our "exactly soluble point" in the 
model. Before we derive the constraints which we encounter 
in tuning the ground state manifold (GSM) to Vd, we must 
first explain Vd and how we map it to dimer coverings. 
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III. d-ISOTOPY AND ITS LOCAL SUBSPACES 



1: 



Although we will eventually be considering a system on a 
lattice, it is useful to begin by defining the Hilbert spaces of 
interest, V d and V d , in the smooth, lattice-free setting. Con- 
sider a compact surface Y and the set S of all multiloops 51 
X cY. If dY (the boundary of Y) is non-empty, we fix once 
and for all a finite set P of points on dY with X n dY = P. We 
assume Y is oriented but X should not be. There is a large vec- 
tor space C 5 , of complex valued functions on S. We say X and 
X' are isotopic (X~X') if one may be gradually deformed into 
the other with, of course, the deformation being the identity 
on dY (see Fig. [3}- 




FIG. 3: Isotopy on an annulus: X ~ X' but X ¥> X". 

We may view the isotopy relation as a family of linear con- 
straints on C s , namely V(X) - ¥(X') = if X ~ X'. The sub- 
space satisfying these linear constraints is now only of count- 
able dimension; it consists of those functions which depend 
only of the isotopy class [X] of X and may be identified with 
C^l, where [S] is the isotopy classes of multiloops with fixed 
boundary conditions. Note that since all isotopes can be made 
by composing small locally-supported ones, so the relations 
we just imposed are "local" in the sense that we will be able 
to implement them with purely local terms in the Hamiltonian. 

Let us go further and define an additional local relation 
which, when added to isotopy, constitutes the "<i-isotopy" re- 
lation. This relation is: 



d x y(x)- x ¥{xuO) = o 



(2) 



It says that if two multiloops are identical except for the 
presence of a small (or, it follows, any contractible) circle, 
then their function values differ by a factor of d, a fixed pos- 
itive real number. In cases of interest to us 1 < d < 2, so 
our function is either neutral to or "likes" small circles. We 
call the subspace obeying all these constraints the aMsotopy 
space of Y (with fixed boundary conditions) and write it as 
V d C C C s . The subspace V d {T 2 ) is still of countable 
dimension, or extensively degenerate, on the torus T 2 . 

It is a remarkable fact (see Refs. l3 ll32l36l) that it is very dif- 
ficult to add any further local relations to d— isotopy without 
killing the vector space entirely. For a a root of unity ^ ± 1 , 
and d — a + a there is such a local relation, but in almost all 
cases the natural inner product on V d fails to be positive defi- 
nite. The physically interesting cases reduce to a = e^l\ k+2 \ 
k = 1,2,3, .. .. We call the corresponding d's "special". 

In these cases the local relations are essentially the Jones- 
Wenzl idempotents: 
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(3) 



k = 2: 
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d 2 -l 
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k = 3: 
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d 4 -3d 2 + 2 
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d 4 -3d 2 + 2 
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1 




d 4 -3d 2 + 2 





= 0, (5) 



see Ref. 30 for a recursive formula. These relations define a 
finite dim ensiona l Hilbert space V d (Y) C V d {Y) C C [ S] C C s . 

In Refs. l3ll32l it is explained that these V d (Y) are the Hilbert 
space for DK^ mentioned earlier, d = 2cos7t/ (k + 2). It is ar- 
gued that a Hamiltonian with GSM corresponding to V d may 
be unstable and collapse under perturbation (for k = 1 or 2 
and under a larger deformation for k > 3) to V d . Very briefly, 
TQFTs such as TiK^, can always be defined as the joint null 
space of commuting local projectors 52 , implying the existence 
of a local Hamiltonian with a spectral gap in the thermody- 
namic limit. Once a Hamiltonian H d has imposed d— isotopy, 
i.e. GSM(H d ) = V d , an extensive degeneracy has been cre- 
ated; the only local way of lifting this extensive degeneracy 
(to a finite degeneracy) without creating frustration 5 ^ is to 
impose the Jones-Wenzl projector as a constraint. Although 
frustration may arise in a Hamiltonian describing a topologi- 
cal phase, we know that these phases can be produced by an 
unfrustrated Hamiltonian. Thus it is an attractive ansatz that 
near H d will be some H d e with GSM(H tx i) = V d . 

To this point our discussion has contemplated smooth mul- 
tiloops X on a surface Y; now it is appropriate to move to a lat- 
tice setting. It is an old idea (see e.g. Ref. l37h to turn a dimer- 
ization (perfect matching) J7 into a multiloop 3^ U J by using 
a background dimerization to form a 'transition graph' . If 
the lattice L is bipartite then %^ U J can be oriented in a nat- 
ural way, leading to conserved quantities 54 - integral winding 
numbers - which are disrespected by the Jones-Wenzl (JW) 
relations. For this reason we pass over the square and hexag- 
onal lattices as poor candidates for the imposition of the JW 
relations crucial to the passage from V d to V d . Thus we con- 
sider perfect matchings on the triangular lattice T . By fixing 
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%, as in Fig.|2j without small rhombi with two opposite sides 
red, as the preferred background dimerization we obtain the 
fewest equations in defining V c / and also achieve ergodicity 55 
under a small set of moves. Unlike in the usual case, the the 
background dimerization, is not merely a guide for the 
eyes, but will be physically distinguished: the chemical po- 
tentials and tunnelling amplitudes will be different for bonds 
of different color. 

Let us list here the elementary dimer moves that preserve 
the proper dimer covering condition: 

• Plaquette (rhombus) flip - this is a two-dimer move 
around a rhombus made of two lattice triangles. De- 
pending on whether a "red" bond forms a side of such 
a rhombus, its diagonal, or is not found there at all, the 
plaquettes are referred to, respectively, as type 1 (or 1'), 
2, or 3 (see Fig.|4). The distinction between plaquettes 
of type 1 and 1' is purely directional: diagonal bonds 
in plaquettes of type 1 are horizontal, for type 1 ' they 
are not. This distinction is necessary since our Hamilto- 
nian breaks the rotational symmetry of a triangular (or 
Kagome) lattice. 

• Triangle move - this is a three-dimer move around a 
triangle made of four elementary triangles. One such 
"fiippable" triangle is labelled 4 in Fig. [4] 

• Bow tie move - this is a four-dimer move around a "bow 
tie" made of six elementary triangles. One such "fiip- 
pable" bow tie is labelled 5 in Fig.@] 

To make each of the above moves possible, the actual dimers 
and unoccupied bonds should alternate around a correspond- 
ing shape. Notice that for both triangle and bow tie moves we 
chose to depict the cases when the maximal possible number 
of "red" bonds participate in their making (2 and 4 respec- 
tively). 




FIG. 4: Overlap of a dimer covering of T (shown in thick black) with 
the red covering corresponding to the special sublattice . Shaded 
plaquettes correspond to various dimer moves described in the text. 
Green sublattice is not shown. 

Note that there are no alternating red/black rings of fewer 
than 8 lattice bonds (occupied by at most 4 non-colliding 
dimers). Ring moves only occur when red and black dimers 
alternate; the triangle labelled 4 in Fig.|4]does not have a Ring 
term associated with it, but the bow tie labelled 5 does. 



Here is the correspondence between the previous smooth 
discussion and rhombus flips relating dimerizations J of 1 . 
The surface Y is now a planar domain or possibly a torus (pe- 
riodic boundary conditions). A multiloop X in Y becomes 
X = %, U 3 (with the convention that the dimers of %, PI 3 be 
consider as length 2 loops or bigons). What about isotopy? 
Move (2) certainly is an isotopy from U 3 to %^ U 3 ' but 
move (2), by itself, does almost nothing. It is not possible to 
build up long motions from (2) alone. So it is a peculiarity 
of the rhombus flips that we have no good analog of isotopy 
alone but instead go directly to d— isotopy. Move (1) must be 
considered in two different forms (1) and (1') as a result of 
the color differences (see Fig. |4). The reader may expect us 
to impose the relations associated with type (5) and (1), (1') 
moves: 

d 3 * [ Jr^j ~ T (\'\) = "• (6a) 

<i xv (/7) l|J (iI") ( 6b ) 

since we pass from one to four loops in do*at and zero to one 
loop in \6b\ . 

However, by choosing, instead a less obvious mapping of 
Vd to functions on {j } we end up imposing one fewer equa- 
tion on the Hubbard parameters. So instead we impose: 




" l| ' (/IT) T (^) ( ?b ) 



and require that a 4 / b = d. 

In other words, the two processes, (1) and (1') for making 
length 2 loops and the one process (5) for fusing four length 2 
loops in one loop of length 8 (see Fig.|4) will have to be tuned 
to satisfy Eqs. @ and a 4 /b = d. 

Given VP e V d , let SP p (j) = a #v F((X Uj7)~) where # is the 
number of length 2 loops of U3 (i.e. #of dimers common to 
3i and 3) and Ujz)~ is the multiloop HtU3 \2i (13, i.e. all 
of the loops except the length 2 loops. This pulls SP, a function 
on multiloops back to a function x Vp(3) on dimerizations. If 
3 and j' are linked by finitely many applications of Eqs. Q 
and the isotopy move, it is now easy to check that *Pp (3 ) = 
a ( # - *)rf( c ' -c )VPp(j / ) where c(resp. c') is the number of non- 
essential loops of length exceeding two in 3^ U 3 ( resp. ^, U 
3'). 

This separate accounting for length 2 loops and the longer 
circles may appear to be a slight of hand. Is there a price 
to pay? In a sense, yes: It is now crucial that a combina- 
torial relation which mimics the smooth JW relation by re- 
lating (3i,J 2 ) for k= 1, (Ji,...,j 5 ) for level 2, (Ji,...j7i 4 ) 
for level 3, etc. . . must not in any of its terms change the 
number of length 2 loops. But this is just an additional (and 
readily achieved) condition on the combinatorial form of the 
JW projector and does not influence the algebraic conclusion: 
that (only) for special d is there a unfrustrated local reduction 
Vd C Vd which could be a stable phase and when d is special 
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Vj is unique. Thus the combinatorial analog of V ( / is functions 
on {j7 } obeying Eqs. @ with a 4 /b = d. 



IV. THE ANALYSIS OF A AND A 

In this section we derive the effective Hamiltonian H : © — > 
£> on the span of dimerizations. The derivation is perturba- 
tive to the second order in £ where e = t' hb /U — tg b /U. Ad- 
ditionally, fy/U = C()£ where cq is a positive constant, while 



o(e) and can be neglected in the second-order calcula- 
tions. (In the absence of a magnetic field all t 's can be made 
real and hence symmetric with respect to their lower indices. 
Also, from now on we set U = 1 for notational convenience.) 

As explained in Section [HI] we then shall tune H to the 
"small loop" value a and the "bow tie" value b with a 4 j b = d. 

We account for all second-order processes, i.e. those pro- 
cesses that take us out of 2? and then back to D (see Ap- 
pendix [B] for technical details). As mentioned earlier, these 
amount to off-diagonal (hopping) processes - "plaquette" 
flips" or "rhombus moves" - as well as diagonal ones (po- 
tential energy) in which a dimer pivots out and then back into 
its original position. The latter processes lead to renormalisa- 
tion of the bare onsite potentials ft, which we have adjusted 
so that all renormalised potentials /},■ are equal up to correc- 
tions o(e 3 ). The non-constant part of the effective Hamilto- 
nian comes from the former processes and can be written in 
the form: 



H 



E(4 



<g>l) 8/ 



(8) 



where H Iy is a 2 x 2 matrix corresponding to a dimer move in 
the two-dimensional basis of dimer configurations connected 
by this move. 8 I} = 1 if the dimerizations I , y 6 © are con- 
nected by an allowed move, 5 Iy = otherwise. 

Therefore it suffices to specify these 2x2 matrices Hj 3 for 
the off-diagonal processes. For moves of types (l)-(3), these 
matrices are given below. Type (1), e.g. rhombus (VY) in 

Fig m 



V 

-2t b t b 



~7t b t b 



V 



rb 



-2c e 2 



-2c £< 



v rb 



a — 1 
-1 a' 1 



(9) 



Type (!'), e.g. rhombus (VW) in Fig. [5] 



"bb 
-?t b t b 



-2t b t b 




v bb 

-2c e 2 



-2c £' 



" rb 



a — 1 
-1 a' 1 



(10) 



Type (2), e.g. rhombus (UV) in Fig.|5] 



"bh 



-Wbb) 



-2(t r bbY 



' bb 



l bb 

-2e 2 



-2e< 



"hb 



1 -1 
-1 1 



(ID 



since isotopic multiloops should have equal weights. 
Finally, type (3), e.g. rhombus (WX) in Fig. |5] 



'bb 



-2{t?-V 



-2(4) 



"bb 



bb) - 



bb 






V bb 








(12) 



provided k > 1, since it represents a "surgery" similar to 
Eq. l|3}- (For k = 1, on the other hand, this matrix must 
be proportional to ( _j "j 1 ) .) As follows from Eq. (0, at level 
k = 1 multiloops which differ by a surgery should have equal 
coefficients in any ground state vector *P while at levels k > 1 
no such relation should be imposed. We use the symbol "~" 
to denote proportional via a positive factor. 

The matrix relations d9H 21 yield equations in the model 
parameters (we consider now only the case k > 1): 



Types (1)&(1') : 

Type (2) : 
Type (3) : 



v gb 



and 



"bb 



2aco£ 2 
= 2a 
= 2e 2 



v bb - 
--v b 

rg 



"bb 



= 



(13a) 
(13b) 
(13c) 
(13d) 



We suppose that the Hamiltonian has a bare ring exchange 

term, Ring in Eq. Q: Ring = ( x 2 ~ C3£ ) for some con- 

\-c3E y J 

stants C3, x, y > 0, and consider the additional equations 
which come from considering Ring as a fluctuation between 
one loop of length 8 (labelled (5) in Fig.|4ji and four loops of 

length 2 . It follows from Eq. ( fTal i that Ring ' '' 1 

so: 



-1 l/b 



,b>0 



a 4 /b — d, x/e 2 = bcs, y/e 2 



(14) 



are the additional equations (beyond Eqs. H31 ) to place model 
A at the soluble point Vj. The justification of the diagonal en- 
tries in the ring term is that the 4 particles in ^0 and ^1 form 
a square which has some cost over an ideal Winger crystal of 
particles. This cost can be influenced by the local chemical 
environment so the entries do not need to be equal, though 
x — y is the most natural case. 

Computer studies 38 show that with the Ring term present all 
J are connected by repeated application of the Hamiltonian 
(and respect the d— isotopy relation). When the Ring term 
is removed, it still appears that any two dimerizations J and 
J7 ' which determine isotopic loops nesting patterns (2?. U J ) 
and (2t Uj') (use our conversion to consider bigons as loops) 
communicate. So, e.g., in a disk if H_ Uj7 and Uj' have 
the same number of loops and the same combinatorial nesting 
relations, and if J , and y' all agree near the disk boundary, 
then moves of type (1), (1'), (2), and their reversals will con- 
nect J and y' '. Since the Jones- Wenzl relations do violence to 
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the nesting patterns it is hard to imagine any additional con- 
served quantities which could allow the ground state without 
the Ring term to collapse to a richer (move degenerate) topo- 
logical phase than <DKk itself (for some k). So we suspect that 
the Ring term is redundant if the goal is to arrive in a topolog- 
ical phase. 

However, there is a caution which should be issued if the 
ring term is omitted. Because of the possibility of defining a 
separate circle value a for bigons and a value b for breaking 8- 

4 

gons, with only the relation L \ = d tying the model to a fixed 
level^the A~ model (Eq. Q with Ring omitted) could still be 
tuned to any level k because only a and not b has been given. 
Thus, for example T)K\ can still arise if an effective Ring term 
with /? = a 4 emerges. 

But, as we have remarked, physically it is most natural to 
assume x — y in the Ring term. This forces b = 1 and a = d x l A . 
Note that we do not want fluctuations on the alternating green- 
black 8— bond rings as this would mix distinct topological sec- 
tors and no such terms are in the Hamiltonian Q. 

Our simplest bosonic candidate for a "universal quantum 

computer" would be model A tuned to a = ^ ; x — 

y = C3£ 2 in Ring (and no green-black Ring term). 



V. MODEL B (SPINLESS BOSON/NO RING TERM) 

For this model we work entirely within the extended Hub- 
bard model, H, as in Eq. Q with no Ring term. What replaces 
the 4-particle ring exchange term is a flip of a new rhombus 
of type [r,r): /•./ (e.g. rhombus QP in Fig.|5jl which can be 
interpreted as an alternating ring of length 4. We create an ex- 
tensive system of (for example, bi-periodic) "defects" in the 
red sublattice . A defect is made by rotating four red edges 
by a 1 /8-turn around a bow tie. Inspecting this defect reveals 
two (r, r) -rhombi adjacent to the defect. The required (r, r)- 
rhombus flip can easily be coaxed out of the Hubbard model 
at o(e 2 ) so there is no need to include an ad hoc Ring term. 
(Without a defect a Ring term does not arise until order O (e 4 ) 
- see AppendixlBl) 

We can construct a d— isotopy space V,/ as the GSM of 
model B in a similar fashion to model A~, but the price of 
having introduced the defect is that there are now many more 
exceptional cases for type (1), (2), and (3) moves. A pro- 
liferation of colors must be defined so that raw potentials /j,- 
when renormalized by hopping in all the different local envi- 
ronments come out equal up to O (e 2 ). We introduce the colors 
(see Fig.|5j: g(green), fo(black), /(indigo), r(red), /(lavender), 
and G(thick green), B(thick black) which are identical in tun- 
nelling properties to g and b (resp.) but require different on- 
site potentials /jq Vg and /jb ^ f*b- 

We take the following relations (which are not the most 




FIG. 5: Sublattice defect, a single bow tie rotation of 3?., in bosonic 
model B. Two rhombi with alternating red and black edges (type 
(r, r)) are shaded lavender. Triangles forming plaquettes whose color 
environment is affected by the defect are labelled A — T and A — T 
(with notations reflecting the obvious inversion symmetry). Trian- 
gles labelled U — Y form regular plaquettes. 

general possible): 

= 4=^=^=4 d5a) 

E = t'bb = tgb = t'bl = t'ii = { bi = { gg = { gG (15b) 

c e = t\ r = tf r = t x hr = t' bb , for x ^ / (15c) 
cie = 4 (15d) 

(Note: With regard to tunnelling amplitudes, / behaves like b 
except t'u = E which has no ^-analog.) 

Near the defect in many special rhombi occur and must 
be considered. All labelled rhombi (those bearing two let- 
ters in Fig. |5J can be classified by the topological effect of a 
two-dimer move: types (1) and (T) create or absorb a loop 
of length 2, type (2) is an isotopy move, type (3) is a surgery 
and an additional type (r, r) (e.g. rhombus {QP)) breaks a 4- 
gon into two bigons (or the reverse). Regard the interaction 
potentials v as variables and for each labelled rhombus write 
the matrix relations, extending Eqs. (f^i — J12i . we may now 
solve for the v's as functions of the t's as shown below J 1 61 - 
I17> . To shorten the list of v's we assume that capitalisation of 
indices for v has no effect. Similarly, we assume that / and b 
are interchangeable. 

Type (3): = vf fe = vf g (16a) 

Type (2) : 2e 2 = vfc = v r bl (16b) 

Types (1)&(1') : 2ac Q E 2 = v b gb = v b bb = v* = v% (16c) 

and 2a- 1 c e 2 = v b rb = v b rg (16d) 

The only relation not previously appearing in model A~ 
comes from from type (r, r) rhombi (e.g. (QP)): 

l-2(4) 2 v[ r ) - {-I b-i) 

=> v' bb = 2bc\z 2 , = 2b- x c\z 2 , (17) 
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with the additional constraint a 2 Jb = d. 

In principle, all labelled rhombi give relations on v's as 
functions of the f's, but we have set things up so that these 
further relations are redundant with the ones coming from VY, 
VW, UV, WX and QP which we have already used in deriving 
dl6!17i . In the following table we begin a verification which 
the reader may complete. 

Type (!'), rhombus (//): 



_2 t b t b^ 



-2t" t" 



v b bb 



(18) 



Type (2), rhombus (I A): 



-2(4) 2 



v bbi 



v h t r 
v bb> l bb 



(19) 



Type (!'), rhombus (AB): 



-2t B t B 

Ll gg l rb 



-2t B t B 

Ll gg l rb 



rg 

t B - 
hb ~ 



c e, C 



(20) 



VI. FERMIONIC MODELS C AND D 



For C the Hamiltonian is again Q but with c, and c\ spin- 
less fermionic annihilation and creation and operators. The 
lattice will be triangular T . The special background dimer- 
ization will be the red sublattice The green sublattice Q 
is also marked and edges of % := T \ (X U Q ) will be called 
black, see Fig.|2] 

The main difference from model A is that we introduce "in- 
direct" hopping between adjacent edges which are angle in- 
sensitive. We regard the vertices of 1 as additional virtual 
sites at a higher chemical potential through which a particle 
can hop. The benefit is that the type 3 move becomes sym- 
metric: the 60°, 60° hops cancel the 120°, 120° hops. We also 
must include some "direct", 60° only, hopping terms to avoid 
killing the desired processes 1 and 2. Here are the tunnelling 
amplitudes (again setting U = 1). 

indirect: There will be a tunnelling amplitude t between any 
two edges of T which share a vertex (regardless of color or 
angle) t = e. 



cie, t b b = c 2 e, t r bb = c 3 £. 



direct: 4 = 

Ring term: ^ 

Here are the terms in the Hamiltonian H and the corre- 
sponding equations: 



Type(l): 



Typed'): 



v gb 
-2cic 2 e 2 



-2cic 2 e 2 



-2c 1 c 2 e" 
v b 



-2cic 2 e" 



"rb 



a — 1 
-1 a' 1 



a — 1 
-1 a' 1 



(21) 



(22) 



Type (2): 



-2c 2 2 e 2 - 



'i,h 



(dirdir) (dir ind) ( ind ind ) 
-4c 2 £ 2 -0 



- 4c 2 e 2 



-2c 2 8 2 



1 -1 
-1 1 



Type (3): 



Type (5): 





4) 


-Co 


S) 











Ring- 



bc4E 2 —C4E 2 



-C4E 



b- l c A e 



(23) 



(24) 



(25) 



We must therefore impose the following equations: 

a 4 /b = d 



(2) : 

(1) and (1') : 

(3) : 



v bb 

„b . 



(2c: 



"bb 



-4c 2 )e 2 
2acic 2 e 2 



(26) 



and v b h = v b g = 2a 'cic 2 e 2 



'hb 



rb 

= 



As in bosonic models, the required potentials can be solved 
for from the tunnelling amplitudes. 

As with model A, we may simply omit the Ring term from 
model C to get model C . Under the assumption that an ef- 
fective ring term with b = 1 emerges, C could be tuned to 
yield •DKk. 

Finally we treat model D. We use the same 3^, red sublat- 
tice with defects as in model B. The necessary colorings are 
illustrated in Fig. [6] 




FIG. 6: Sublattice defect in fermionic model D. 

We list the tunnelling amplitudes: 

indirect: t = e between all adjacent edges, e,- n ej ^ 0. 

direct: 



= 


t 8 - 
l bb — 


t X 


Cl£ 


-t b 
— l rb 


-t b 

— l lb 


C 2 £ 


-t b 

- tgb 


= t h 
gg 


C36 


— hb 


= tlb 


C4E 


-J 

— hb 





(27a) 
(27b) 
(27c) 
(27d) 
(27e) 
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We have only listed colors b,r,g,l but some of these may ap- 
pear in two or three forms because of differing local envi- 
ronment. Though each of these forms needs an individual fj 
(chemical potential), the t and v symbols will depend only on 
the 4 given colors. 

Near the defect, to maintain the equations of model A we 



need to set: v 8 , — vf t , v r — • ■' • '" 



v 'bb- v " g i = v i;,< V rl = V rb> A = v bb- 



gb ~ v bfr v bl 

This effectively maintains the constraints associated to 
rhombi of types 1,1,2, and 3 near the defect. 

The lavender rhombi (Fig. [6) yield a new term which must 
be constrained by the equations below. 



v bb 
-2c 2 e 2 



Me 2 



b -1 

-1 b- 1 



(28) 



so we have the equations: 



v bb = Ibciz 



2„2 „/ 



4 C ' 



2b- l clt 2 , — =d 1 (29) 
b 



in addition to the equation from model C: 

v r bl = v b bb = (2c 2 2 + 4c 2 )e 2 
v b g i = v b gb = v b bb =2ac lC2 e 2 



v° rl =v D rb = 2a- l cic 2 e 



(30a) 
(30b) 
(30c) 
(30d) 



To find solutions it is only necessary to choose a and b com- 
patible, pick the positive constants c\, ■ ■ ■ ,ca, and compute the 
potentials v. 
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APPENDIX A: DEFINITION AND PROPERTIES OF <DK k 

The information in this appendix allows one to predict the 
outcome of any Aharonov-Bohm experiment conducted on a 
material thought to be in the topological phase DK^. 

There is a well known 39 family of topological quantum 
field theories (TQFTs) SU^)*. These are the unitary SU(2) 
TQFTs at level k. In the quantum group approach 40 in order to 
define precisely this theory (and not e.g. its time reverse) one 
sets the deformation parameter A for the universal enveloping 
algebra of su(2) to A = e 2lti /4(i+2). 



SU(2)jt has has particles species indexed by SU(2) repre- 
sentations (or 'isospins') j = 0,1/2, l,...,k/2. These k+ 1 
irreducible representations are subject to the following rules: 

The fusion space Vjy = C if 



0<ji,ji,j3,<k/2 

71 +72 > 73, 
73 + 71 >72, 

72 + 73 >7i> 

7i + 72 + 73 < k. 



(Ala) 
(Alb) 
(Ale) 
(Aid) 
(Ale) 



= otherwise. 
The S— matrix is 



Sjoz = sin ( 2 >i + 1 ) ( 2 >2 + 1) */(*+2)) , 

0<i,j<k/2. (A2) 

Th ere is a variant of SU(2)i worked out in 41 (also see 
Refs. l32l3rt) based on the Kauffman bracket for links but now 
with A = ie 2,ti /4(A:+2) in the relation of Fig.0 



A 



A 



FIG. 7: Kauffman bracket 

On a torus T 2 , K^T 2 ) would be a subspace of the functions 
on links in the solid torus, respecting the above relation, the 
level k Jones- Wenzl relation, small unknots = d — 2 cos n/k + 
2, and finally the relation that a right handed kink is equal 
— A~ 3 . Kk also has k + 1 orthogonal irreps 0, . . . ,k with the 
same fusion relations as SU(2)jt but the S— matrix is Suh — 

(-l) 4hh ^(M2ji +l)(2j 2 + l)7l/(fe + 2)). 

Kk is a Chern-Simons-Higgs theory in which two SU(2) 
Chern-Simons gauge fields, with levels k and 1, are tied to- 
gether by a condensate which transforms in the isospin k/2 
representation of the first and the isospin 1 /2 representation 
of the second: 



S =fc5cs[fli]-5cs[fl2]+5Hig gs [*,fli,fl2] (A3) 



where 



■s'cs — / e^P ( a^ v ap + \ fabcafava? 



1 

4jt 



1 r f ,2 

— / tr aArfflH — flAflAa 
471 J V 3 



(A4) 



and 



5Kgg.[*,fli,fl2] = / | (i3„ + a%Tfj +4 M 4) fjjf+Vi®) 

(A5) 

a = 1,2,3 are su(2) Lie algebra indices; T-, t- are su{2) gen- 
erators in the isospin k/2 and 1/2 representations, respec- 
tively. V(<S>) is minimized at non-zero |<t>|. The presence of 
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the condensate restricts the possible particle types so that half- 
integer isospins under SXJ(2) k are necessarily spin-1 /2 under 
SU(2)i. This leads to the additional minus signs in the braid- 
ing statistics and the 5— matrix. 

For k even, det 5=1 but for k odd det 5 = 0. This singular- 
ity is explained by the isomorphism: 



1 V 

.1 1, 



5 integer, integer 



for k odd . 



(A6) 



Thus K k is only a true TQFT for k even. Physically, the 
problem with odd k is that it is not possible to distinguish the 
isospin-j particle from the isospin-(£/2 — j) particle by any 
braiding operation 3 ?^. 

Planar arc diagrams lead both to K k , via Fig. [7] and to 
SU(2)fc, via the Rumer- Teller- Weyl theorem on SU(2) rep- 
resentations (see Ref. 42 for an exposition of both). The sub- 
tle difference is that certain — 1 's, Frobenius-Schur indicators, 
occur in the representation theory coming from the fact that 
the fundamental representation of SU(2) is antisymmetrically 
self-dual, whereas these signs are absent in Kauffman's the- 
ory. 

We come now to T>K k , the Drinfeld double of K k . It is 
a subspace of functions on links in surface x interval (same 
relations as listed for solid torus) and has its own naturally 
defined positive definite hermitian inner product and a unitary 

5— matrix, 5 41,43 . When k is even 

'DK k =Kl(g, c K k <is Hilbert spaces, and 5 = 5 <g) 5, (A7) 
when k is odd, we still have: 

VK k =K* k ® c K k , 

but now 5 = T <X> 5i nte ger, integer ® ^integer, integer: (A8) 

i.e. 5 is the tensor square of the nonsingular even part of 5 
tensor T, the 4x45 matrix of the level k — 1 theory T>K\, In 
Kitaev's "toric code" notation 28 the irreps of T>K\ are: 0,e,m, 
and em; the trivial, electric, magnetic, and electric and mag- 
netic excitations respectively. The Lagrangian for this sector 
is that of Z2 theory 4 ^. 
In this basis: 











(A9) 



m 
em 



'1/2 1/2 1/2 1/2 

1/2 1/2 -1/2 -1/2 

1/2 -1/2 1/2 -1/2 

J/2 -1/2 -1/2 1/2 

The isomorphisms in Eqs. JA6> - (IA8> can be made explicit 
as follows (see Refs. l32l36l) . 

K k has basis 0, . . . ,k/2. Define i — k — i, < i < K. If k is 
odd then exactly one of i and i is even; denote that one by i e . 
So when k is odd we write i = (i e ,Oi) where a, = + or — as i 
is even or odd so (i e ,+) = i e , (ie,—) = ie- 

T>K* k (g) DK k has basis (i,j,l,m), < i,j,l,m < k/2, and 
(i,j,l,m) = (i e ,Gi,j e ,Gj,l e ,ai,m e ,a m ) in the case k = odd. 
Declare an isomorphism 9 from (a,, a,) and from (o;,O m ) 
to the toric code excitations by: (+,+) — > 0,(— ,+) — > 
e, (+, — ) — > m, and (— , — ) — * em. Then to realize ( IA8i . map 
(i,j,l,m) to (Q(G i ,Gj),Q(Gi,G m ),i e J e ,le,m e ). 



APPENDIX B: NEARLY DEGENERATE PERTURBATION 
THEORY 

Although for the purposes of this paper we do not use 
anything beyond the second-order perturbation theory, here 
we present a general scheme useful for extending this type 
of analysis to higher orders. In particular, three and four- 
dimer moves will appear as higher order terms in the effective 
Hamiltonian. 

The perturbative scheme presented here closely follows 
that developed in Ref. |45j The idea is to recursively block- 
diagonalise the Hamiltonian, order by order eliminating terms 
that change the number of dimer collisions (a collision oc- 
curs when two dimers share a vertex). The resulting effec- 
tive Hamiltonian then contains terms (up to a given order) that 
only connect states within sectors with a fixed number of such 
collisions. To proceed with this programme, we rewrite the 
original Hamiltonian Q as 



H = T + 7i + + T x + T-x +D + W. 



(Bl) 



Here T +Ti + T- 1 +T x +T-x = - Uj [c\ cj + h.c.),D = 

YtiHifli +D(/J)eix,£0 Vij n i n j is me "low-energy" part of poten- 
tial energy, while W combines the remaining "high-energy" 
terms in H, i.e. all dimer collision interactions. To, T\ and 7Li 
represent the dimer moves that respectively do not change, in- 
crease by one, or decrease by one the number of collisions. 
T x results in a dimer ending on top of another dimer (dou- 
ble occupancy in the original particle language), while T—x 
undoes that. Since we have already chosen Uo — °° thus re- 
stricting our attention to the single-occupancy subspace within 
which T±x have zero matrix elements, we should drop these 
terms from our consideration. Notice that T^ = T- m and 
[T m , W] = -mUT m , m = 0,±1. 

Our strategy is to recursively construct the operator i5 such 
that a "rotated" Hamiltonian H given by 

H = e is He- lS = H+ [iS,H] + — [iS, [iS,H]] + ... (B2) 

is block-diagonal as described above. As a 

first approximation, we choose i5^ = (7i — 
T-i)/U which leads to H^ 2) = T Q + D + W + 
{[T h D}-[T- U D} + lT h To]-[T- h To\-[T-i,Ti})/U + 
(x 3 /t/ 2 ), where x = max{fy,Vjy,///} (here we follow the 
numbering convention of Ref. 45 where ffW = H). This 
procedure eliminates T\ and 7Li but generates a slew of 
smaller terms, all of which, except —[T-i,Ti]/U, change the 
number of dimer collisions. The next step is to correct i5 in 
order to eliminate these new unwanted terms, repeating this 
procedure recursively up to any given order. The number of 
terms rapidly gets out of hand, and we used Mathematica to 
keep track of them up to the fourth order (i.e. keeping terms 
of order x 4 /t/ 3 ). 

We are particularly interested in the terms that connect 
states within the lowest energy, zero collisions sector. Thus 
we can write H = H' + H" where H" consists of terms van- 
ishing in this sector (i.e. the matrix element (oc|Z/"|p) = for 
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any two proper dimer coverings |a) and |p)). Then the effec- 
tive low-energy Hamiltonian is given to the fourth order by: 

ffW = D-T- 1 T 1 /U 

X - DT i Ti + T \DT\ - i 7Li7iD + r_ 1 7b7i J /U 2 



1 



■DDT^iTi +OT i£>7i - T-\DDT\ 



+ T-\DT\D— — T-\TiDD + DT-iTqT\ 
- T iDToTi - T i TqDTi + T i TqT\D 

T^TqTqT! + TLiT-iTLiIi - l - T^T^Ti] /U 3 



D 



TLiTi 1 



r_ir ri 



_{?;_, [D,r 1 ]+h.c.} + -^ 



[D, [D,7i]] +h.c.} — {T-\Tq [D,T\] +h.c.} 

- t i Tbror! + r_irir_iri - i r_ir_irir^ /t/ 3 (B3) 

The advantage of rewriting this Hamiltonian using com- 
mutators such as [D,7i] becomes clear if we recall that D = 
Eifiifii + £((j)£x,^0 Vij n i n j i s diagonal in the local position 



basis. This allows us to combine the terms as follows: 



0G 



■r_iri + -(r_i[A7i]+h.c.) 
-i(r_! [d,[d, ri]]+h.c.) 



1 (y ( d n -d§ 



P) 



d n — d^ 
U 



x(a|r_i|«)<n|ri|P) + (a«->pU (B4) 



and 



(a | (T-iToTi - (T^iTo [D, Ti] +h.c.))| P) 



U 



-(a<-p) (B5) 



where \n) is an excited state with one dimer collision and d n , 
dp are the eigenvalues of D. Written in form JB4HB51 . these 
terms explicitly depend only on local differences in the two 
dimer configurations rather than the expectation values of in- 
finite sums D. 
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